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Abstract 

We derive the pion mass and the elastic pion-pion scattering amplitude in the QCD 6- 
vacuum up to next-to-leading order in chiral perturbation theory. Using the modified inverse 
amplitude method, we study the ^-dependence of the mass and width of the light scalar meson 
a( 500) and the vector meson p(770). 


1 Introduction 

Light flavor quantum chromodynamics (QCD) with up, down and strange quarks is a fascinating 
theory, that features only a few parameters. Variations of these parameters have been explored 
in great detail in the last decade. We mention in particular varying the quark masses to make 
contact with lattice simulations, but also variants of the theory with more light quark flavors 
or different numbers of colors have been studied. These studies mostly focused on the fate and 
exploration of the dynamical and explicit symmetry breaking that light flavor QCD exhibits A 
much less explored territory relates to the 0-term of QCD, 

(1) 

where a = 1, ..., 8 are color indices and G ^ is the gluon field strength tensor. While the upper 
limit on the neutron electric dipole moment poses a stringent limit on the value on 0, where the 
latest determination gives |0| < 7.6 • 10” 11 [1], it has been argued that based on the string theory 
landscape, 0 might take natural values and that it is difficult to achieve a tiny value for it, see e.g. 
Ref. [2]. Independently of that, it is interesting to explore QCD at larger values of 9. QCD at 9 ~ n 
was already investigated in Refs. [3, 4]. Also, the pion and the nucleon mass were calculated at 
fixed topology in Ref. [5] to explore the connection between lattice QCD and physical observables. 
Further, Ubaldi [6] studied the effects that a non-zero strong-CP-violating parameter would have 
on the deuteron and di-proton binding energies and on the triple-alpha process using a somewhat 
simplified nuclear modeling. Even so the relevant energy scales in these systems exhibit some 
fine-tuning, no dramatic effect of varying 9 was found. We also note a recent study on the relation 
between confinement and the 9- vacuum, see Ref. [7]. 

The 0-dependence of the pion mass has been given at leading order (LO) in chiral perturbation 
theory (ChPT) [5]. Here, we want to derive the pion mass in the 0-vacuum up to next-to-leading 
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order (NLO) as well as the corresponding pion-pion (7T7 t) elastic scattering amplitude. Further¬ 
more, we will calculate the 0-dependence of the lightest resonances in QCD, the scalar meson 
ct( 500) and the vector meson p(770). These are not only interesting by themselves, but also are 
important components in precision modelings of the nuclear forces, as used e.g. in the work on ex¬ 
tracting limits on variations of the Higgs vacuum expectation value from the element abundances 
in Big Bang nucleosynthesis [ ] (for related works using different frameworks, see Refs. [9, 10]). 
In the following, we assume that the a and the p are generated from a resummation of pion-pion 
interactions evaluated at NLO in ChPT in the 0-vacuum. Having calculated the 0-dependent pion- 
pion scattering amplitude, it is straightforward to implement it into a unitarization (resummation) 
scheme that can generate the light resonances and thus gives their 0-dependent properties. To be 
specific, we make use of the so-called (modified) inverse amplitude method, which is one, but not 
the only available, unitarization scheme that allows to perform this task. 

Our work is organized as follows. In Sec. 2 we discuss the chiral effective Lagrangian in the 0- 
vacuum, with particular emphasis on the two-flavor formulation and also including strong isospin 
breaking. Next, the 0-dependent 7T7 t scattering amplitude is constructed at NLO in Sec. 3. Armed 
with that, we then come to the central section 4, where the mass and the width of the a and the 
p are calculated as a function of 0. We end with a short summary and outlook in Sec. 5. The 
appendix contains some discussion of the vacuum alignment at NLO. 

2 Chiral effective Lagrangian in the 0-vacuum 

At the lowest order, O ( p 2 ), the SXJ(N) chiral Lagrangian in the 0 vacuum is [11, 12] 

C 2 = F t {D^D^U)+ F t ( X U^ +x'U), (2) 

where F is the pion decay constant in the chiral limit and x = 2HA4 exp (iO/N). Here, A4 is the 
real and diagonal quark mass matrix, and the low-energy constant (LEC) B = E/E 2 , with E the 
absolute value of the flavor-averaged quark condensate in the chiral limit. The field U £ SU (N) 
collects the Goldstone bosons of the theory. However, only for 0 = 0 the vacuum expectation 
value of it, which is the solution of the equations of motion for the zero-momentum mode, is 
trivial, i.e. Uq = 1. In the general case, the vacuum is shifted from the unit matrix, and the 
vacuum alignment can be determined by minimizing the potential energy. Therefore, it is useful 
to separate the ground state t/o from the quantum fluctuation U containing the Goldstone boson 
fields as U ( x) = U 0 U(x). Thus, the ground state of the theory is given by minimizing the potential 
energy 

v 2 = -1 ( (c/ 0 f e W/N + U 0 e- ie/JV ) M ) . (3) 

Because A4 is diagonal, Uq can be taken as diagonal as well without loss of generality. The special 
unitary matrix Uq can be parametrized as 

U 0 = diag [e llpi , e^ 2 ,..., e lVN } , ^ <pf = 2nn {n £ Z) . (4) 

/ 

This leads to 

V 2 = -Y,Re(e- ie/N UvM) = -E^cos (tp f (5) 

with rrif the quark mass of flavor /. 

For 9 = 7r, because the theory is periodic in 0 with a period 2-7 t, the Lagrangian is invariant 
under CP and P transformations as they change 0 to —0. However, it is well-known that at 0 = n, 
there is a spontaneous CP breaking, called Dashen’s phenomenon, because the CP conserving 
stationary point of the action is in fact a maximum and there are two degenerate CP violating 
vacua which are obtained by minimizing the potential energy [13]. 
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2.1 Two-flavor case without isospin symmetry 

In the rest of the paper, we will consider the two-flavor case and use the following parametrization 

U 0 = diag{e* e"*}, U = e^ F , = -L . (6) 

We see that the angle ip and the neutral pion field always appear in a linear combination tp + ir 0 /F. 
Therefore, finding the stationary solution for XJq by minimizing the potential energy with respect 
to ip is equivalent to removing the tree-level tadpole for the neutral pion [13]. The minimization 
of Vi gives [5] 

, 9 9 

(m u + rad) sin ip cos - — (m u — cos ip sin - = 0 

Q 

<t=> t&iup = — etan - , (7) 

where the average light quark mass fh = (m u + ma )/2 and the parameter e = (mj — m u )/(2fh), 
that quantifies strong isospin breaking, are introduced. Actually, at 9 = n and e = 0 the Eqs. (7) 
do not depend on ip at all. This leads to a paradoxical situation of continuous vacuum degeneracy 
discussed in [3, ] and resolved in [14] taking into account terms of the NLO chiral Lagrangian. 

In the present work we also consider the NLO chiral Lagrangian [11]. In the SU(2)xSU(2) 
notation of e.g. Ref. [ ] it reads 

U = l j(D^D»U) 2 + l j(D ll U'D„U)(D^D I 'U) + ^(x i U + xrf) 2 

+j {D^D^U + D^U') + | (tfF*UF L ^) 

+ l ~Y (f^ud^ + f f v d^d v u) - ^ (x'u - x u ] ) 2 

+ h 1 + h 1 + h 1 -_h , Re (det%) _ h2 ( F L vF L^ + pRpR^) , (8) 

where 

D.O = d^O - ir^O + iOl „ for 0 = U, X 

F^ = d»r v - d v r» - i[r", r v ], F£ v = d^F - d v F - i[F, F] (9) 

with r M and F the right-handed and left-handed external fields, respectively. 

In principle, with the introduction of the NLO Lagrangian as well as the one-loop contribution, 
the vacuum energy has changed, and the vacuum alignment needs to be re-determined. In par¬ 
ticular, the F term in the NLO Lagrangian is not minimized by the LO solution given in Eq. (7). 
This means that the F term induces a shift to the LO vacuum alignment. However, as shown in 
Appendix A, this shift does not affect the calculation of the pion masses and the irn scattering 
amplitudes up to NLO. Therefore, it is sufficient to consider the LO vacuum alignment in Eq. (7) 
for our purpose 1 . 

2.2 ^-dependence of the pion mass 

Substituting Uq with ip given by Eq. (7) into the LO Lagrangian, we get the LO pion mass squared 
in the 0-vacuum [5] 

M 2 {9) = 2Brhcos ^ 1 + e 2 tan 2 ^ , (10) 

which is the same for the neutral and charged pions. 

1 We will not discuss the complications at 6 = n. In that case, one needs to include the I 7 term as it determines 

the whole dynamics [3]. 
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At NLO, the pion masses receive contributions from both one-loop diagrams and the Z 3 and l 7 
terms. The divergence in the one-loop diagrams cancel exactly with that from Z 3 . We obtain 


m 2 + (o) 

M 2 Ad) 


M 2 (9) 


M 4 (9) 


F 2 


32tt 2 


In 


M 2 (0) 

y 2 


+ 2 l r 3 + 2 1 7 


(1 — e 2 ) tan(0/2) \ ' 
1 + e 2 tan 2 (9/2) ) 


M 2 + (9)-2l 7 


M 4 (0) 


P 2 cos 4 (0/2) (l + e 2 tan 2 (0/2)) 


2 ’ 


(ii) 


where Z 3 is the scale-dependent finite part of Z 3 . At 9 = 0, these expressions reduce to the standard 
SU(2) relations derived in Ref. [ ]. Using the positivity bound for l 7 obtained in Ref. [3], we find 
that the charged pion is always heavier than the neutral one. 

For easy reference, we give the corresponding formulae for the much simpler isospin symmetric 
case with m u = nid = fh. In this case, the stationary solution of the vacuum energy has ip = 0. 
The pion mass up to NLO in the 0-vacuum has one additional term compared with that in the 
0 = 0 case, and is given by 


M 2 (0) = M 2 (0) 


M 4 (0) 

F 2 



In 


M 2 (0) 


+ 2Z 3 + 2 1 7 tan 2 


( 12 ) 


with isospin symmetric LO pion mass 

n 

M 2 (9) = 2Bm cos - . (13) 

One sees that even in the isospin symmetric case, the NLO pion mass depends on l 7 , and this 
additional term vanishes at 0 = 0. 


3 7T7T scattering amplitudes in a ^-vacuum 


The 7T7T scattering amplitude at NLO is the building block to generate the light mesons cr(500) and 
p(770) via unitarization. To be specific, we calculate the amplitude A(s, t , u) = A 7r + 7r -_ >7r o, n .o (s, t, u ) 
which is used to get the following combinations with definite isospin (/ = 0,1, 2) 

T°(s,t) = 3 A(s,t,u) + A(t,u,s) + A(u,s,t), 

T 1 (s,t) = A(t,u,s) — A(u,s,t) , (14) 

T 2 (s,t) = A(t,u,s) + A(u,s,t) ■ 


Later, we will also need the partial-wave projection for given isospin I and angular momentum L 


T[(s) 


1 1 
32^2 


r +1 


dz T I {s,t)P L (z) 


1-1 


(15) 


with Pl{z) the pertinent Legendre polynomials. 

Up to the order 0(p 4 ), there are several contributions to the nn scattering amplitude, as 
shown in Fig. 1. Diagram (a) gives 

A( a )(s, t, u ) = 7 ^ |3s + Mg — 2 [M 2 o(0) + M 2 +(0)] j , (16) 

in terms of Mg = M 2 (0), Af 2 + (0) and M 2 O (0) given in Eqs. (10) and (11), respectively. Dia¬ 
gram (b) gives 


A (b )(s,f,u) = 


32 l 7 B A m A 


h 


(s - 2M 2 ) 


4Mg ( s - 2Mg ) +t 2 + u 2 


+ 


3F 4 


F 4 l 

[(1 - e 2 ) 2 sin 2 0 — 2e 2 ] . 


3 F 4 9 


(17) 
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Figure 1: Feynman diagrams for the 7T7 t scattering amplitude up to O (p 4 )- Here the filled circles 
and square denote the vertices from the LO and NLO Lagrangian, respectively. The t and u 
channel two-point loops are not shown. 


Diagram (c) includes the tadpole vertex correction from both the neutral and charged pions, and 
its contribution is 

( 18 ) 


A (C)(s,t,u) = (31 Mg - 20s) A 0 (M e 2 ) 


Ao(m 2 ) is the one-point loop integral (the tadpole) in d space-time dimensions 


H 0 (m 2 ) := ifj , 4 d f 


d d l 


1 


(27r) d l 2 — m 2 


(19) 


with fj, the scale of dimensional regularization. 

The two-point loops, diagram (d) and the corresponding t- and u-channel crossed diagrams, 
give 

A (d)(s,t,u) = |^9 

1 


+ 


+ 


+ 


+ 


6.F 4 

1 

487T 2 

6 F*' 

1 

487T 2 


9 (M 4 - s 2 ) B 0 (s, M 2 , M 2 ) + 10s (m 2 

2Mq + 2Mq (t - 2 u) + t(u - t)] B 0 (t, Mlti-lf) + |(t - 3u)H 0 (m|) 


(s -u)(t- 6M|) | 


2 M$ + 2Mg(u - 2 1 ) + u(t - u) Bo (it, M 0 2 , M e 2 ) + ?(u - 3 t)A 0 (iVf' 


(s - t) (u - 6 M^j | 


( 20 ) 


where the first line corresponds to the s-channel charged and neutral pion loops, the second and 
third lines correspond to the t-channel loop, and the last two lines are for the u-channel loop. Here 
Bo is the scalar two-point loop integral 


B 0 {q 2 := ifi 4 


d d l 


1 


(27r) d ( l 2 — m\)[(l + q ) 2 — m 2 } 


( 21 ) 


We also need to take into account the wave function renormalization for all external lines which 
is represented by diagram (e). This amounts to 


H (e) (s, t, u) = - (2 5Z n+ + 2 6 Z^) (s - M 2 ) = _H 0 (m^ ) (s - M 2 , , 


( 22 ) 


where SZ T = Z v — 1, with the wave function renormalization constant for both the neutral and 
charged pions given by 


Z-„ — 1 


3F 2J 


(m 2 ) . 


(23) 


Using dimensional regularization for the loop integrals and summing up all contributions, we 
obtain a UV divergence-free and scale-independent amplitude. The I 3 and I 7 terms in Eq. (17) 
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cancel with the same terms in Eq. (16) that enter through the NLO mass expressions for the 
neutral and charged pions. The full amplitude reads 


A(s,t,u) 
B(s, t, u ) 

C(s,t,u) 


s - Af 2 


u ) + C ( s ’ u ) > 

-4r(3(s 2 - Mq)J (s) + [i(i - ti) - 2M| t + AM'l u - 2M?1 J(t) 
6-r 4 l L J 

u{u — t) — 2 Mg u + AMg t — 2 Mg J(u) j , 


(24) 


1 


96t r 2 E 4 


^2 ^ — 2 Mq^J + ^2 e — — ^ [ s “ + (t — uf] — 12 MgS + 15Mg j- . 


Here, the Ug are the scale-independent but quark mass-dependent, and thus ^-dependent, LECs 
which are related to the renormalized ones as 


7i 

32tt 2 


|' hs + In 




The finite loop function J is given by 


A*) 


i 

167T 2 


^u(s) In 


g (g) ^ 1 

cr(s) + 1 



As) 



(25) 


(26) 


We find that the tttt scattering amplitude up to NLO in a 9- vacuum, Eq. (24), takes exactly the 
same form as the well-known one in the vacuum with 9 = 0 [ ], and the only change is to replace 

everywhere M 2 ( 0) = 2 Bfh by M| given by Eq. (10). The reason is that vertices from terms of 
the form (y}U + xPA can always be written in terms of Mg, while the Z 7 term from diagram (b) 
gets cancelled with the one from diagram (a). Such a property does not hold at higher orders. For 
instance, considering the 7 T 7 t scattering at O ( p 6 ), there can be a one-pion exchange diagram with 
two CP-violating three-pion vertices (see Appendix A), which does not have any correspondence 
at 9 = 0. This behaviour of the 7 T 7 t scattering amplitude is reminiscent of the Kaplan-Manohar 
transformation [16], which is an accidental symmetry of the chiral Lagrangian at NLO. 


4 ^-dependence of the a and p in the isospin limit 


The cr(500) and the p(770) are the lightest two-flavor non-Goldstone mesons. They can be obtained 
from the chiral perturbation theory amplitudes by militarization. There are various such lineariza¬ 
tion schemes on the market, like the inverse amplitude method (IAM) to be used here [17]. In most 
cases, such a unitarization procedure amounts to a resummation of a certain class of diagrams to 
ensure exact two-body unitarity, which is only perturbative in ChPT, but such resummations are 
usually at odds with crossing symmetry. We do not want to enter a more detailed discussion on 
these issues here (see e.g. the early work in Ref. [18]), but rather employ the IAM as a tool to 
generate the light mesons from the 0 -dependent pion-pion interaction, which automatically leads 
to 0 -dependent properties of the a and the p. 

The scattering amplitude for a given channel (with fixed isospin and angular momentum) up 
to NLO in the IAM is given by 


T(s) 


(W ) 2 

T m {s)-T {A] {s) ’ 


(27) 


where T( 2 )(s) and T^(s) are the irn scattering amplitudes of leading and next-to-leading chiral 
order. This form is valid in the channel with I = J = 1 pertinent to the p-meson. As pointed 
out e.g. in Ref. [19], it requires modification in the I = J = 0 channel due to the presence of 
Adler zeros in the S- wave. The associated unphysical poles can be cancelled in rather natural 
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Figure 2: The 0-dependence of the masses (left panel) and widths (right panel) of the a (blue, 
full) and the p (red, dashed). 


way as derived in Ref. [20] which is called the modified inverse amplitude method (mlAM). The 
corresponding scattering amplitude reads 


T(s) 

A mlAM (s) = T (4) ( S2 ) - 


cwr 


T( 2 )(s)-T (4 )( S ) + A mIAM ( S ) ’ 

0 2 - S A )(s - s 2 ) ( T' (2) (s 2 ) - r ( ' 4) (s 4 )) 


(28) 


S - SA 


where sa denotes the Adler zero of the full partial wave defined by the condition T(s a ) = 0 , 
and ' denotes differentiation with respect to s. The approximative Adler zeros at LO and NLO 
correspond to the energies s 2 and s 2 + s 4 , determined by T( 2 )(s 2 ) = 0 and T[ 2 )(s 2 + s 4 ) + T( 4 )(s 2 + 
s 4 ) = 0, respectively. This mlAM has been used e.g. in Ref. [2 ] to study the quark mass 
dependence of the sigma and the rho. In particular, we will use the LECs l\ and l 2 as determined 
in that paper at the scale p = 770 MeV, 


l\ = (-3.7 ± 0.2) x 1CT 3 , l r 2 = (5.0 ± 0.4) x 10 “ 3 . 


(29) 


As mentioned in Ref. [21], the results of the IAM is insensitive to the values of 1 3 and / 4 as 
long as they are within the uncertainties: 1 3 = (0.8 ± 3.8) x 10 -3 and l\ = (6.2 ± 5.7) x 10~ 3 . 
Taking the central values and the measured values of the pion mass and decay constant M^ = 
138.04 MeV (isospin averaged) and F n = 92.21 MeV, the standard ChPT one-loop expressions 
yield M = 139.46 MeV and F = 86.43 MeV. 

The masses and widths of the p and <7 resonances can be obtained by searching for the poles 
in the complex s-plane of the unitarized amplitudes. When the resonances are above the 7 T 7 t 
threshold, which is the case for the physical pion mass, we need to search for poles in the second 
Riemann sheet. The corresponding pole positions read 


y/sa = (443.1 — i 217.4) MeV, y/Sp = (751.9 — * 75.4) MeV. (30) 

We note that the mass of the p comes out somewhat below the physical value as it is common in 
such unitarization procedures. For a more detailed discussion on this issue, see e.g. Refs. [22, 23]. 


When isospin breaking is neglected, m u = rrid , the vacuum is not shifted, and we can use the 
usual ChPT Lagrangian and amplitudes directly. All the 0-dependence of physical observables 
enters through Eq. (13) which finally leads to the mass and width of the a and the p as a function 
of 0, shown in Fig. 2. The 0-dependence of the <r mass is stronger than the one of the p mass since 
the former is in a S'-wave while the latter is in a P- wave. Also, both widths show a somewhat 
stronger dependence on 0 which is due to the enlarged phase space as the pion mass decreases 
from its physical value when | 0 | increases from 0 . 
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5 Summary and outlook 

In this paper, we have studied the 0-dependence of the lightest resonances in QCD. For that, 
we have derived the charged and neutral pion masses and the pion-pion scattering amplitude 
at NLO in the 0-vacuum. We found that the NLO contributions proportional to I 3 and I 7 and 
entering via the pion mass formula and 7T7r-contact terms cancel each other exactly. The er and 
the p have been obtained from a unitarization of this amplitude using the so-called (modified) 
inverse amplitude method. This automatically generates 0-dependent masses and widths of these 
resonances. Although the pion mass vanishes at 0 = n at LO, no dramatic effects on the masses and 
widths of the a and the p were found. However, it still remains to be seen how such modifications 
change the properties of nuclei, as the nuclear binding is fine tuned, and thus more sensitive to 
such parameter variations. 
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A Vacuum alignment at NLO 

In this appendix, we will discuss the vacuum alignment at NLO induced by the presence of the 
counterterms (/,-terms) and the one-loop contribution. We will calculate the vacuum alignment 
perturbatively. 

Because the shift of the LO vacuum alignment given in Eq. (7) is caused by the NLO terms in 
the chiral expansion, we assume that the angle ip in U$ = diag{e* v , e _l<|P } can be split into 

tp = (p 0 + atp 1 (A.l) 

with ipQ determined by aligning the vacuum at LO, and aipi is the shift coming from the O (p 4 ) 
contribution to the vacuum energy. Here, a is a chiral scaling factor to make explicit that a p>\ 
is one order higher than <p$ in the chiral expansion. 2 The vacuum energy density up to NLO is 
given by 


evac = ~ {x'Uq + X ul) - (x f t/o + xUl } 2 + ^ (x + E/o - X^)' 

-^ + ^3 ^ _ hi_hs Re (det%) + eWioop) ; (A.2) 

where we have neglected those hi- terms which are independent of Uq, and eii c lo ° P ^ is the 1-loop 
effective potential whose explicit expression is [24] (see also Ref. [25] expanded up to 0 4 ) 


e va~c° op ' ) = 3 M\ 6 ) 


A 1 


2 128t r 2 


1 - 2 In 


M 2 (0)\ 

p 2 J 


(A.3) 


We may decompose the vacuum energy density into the LO and NLO contributions with the 
NLO one including all terms proportional to a 


e = e (2) , (4) 

c vac — c vac W u. e vac . 


(A.4) 


2 The introduction of the scaling factor is only for convenience. It will be set to a = 1 after </?i is calculated. 









/ 


/ 


7T U \ 

- 


" - -0 


7T° \ / 

- 


Figure 3: One-loop diagrams for producing a neutral pion from the vacuum. Here (g> and the 
black dot denote CP violating and conserving vertices from the LO Lagrangian, respectively. 


Substituting Uq = diag{e l¥3 ,e lv } and Eq. (A.l) into Eq. (A.2), we get 

4 ac = M^ = 2Bfh^cos^cos(p 0 -esm^sm(po\ , 


e (4) = 

^vac 




sin - cos po + e cos - sin p 0 ) , (A.5) 


where we have neglected the terms independent of p. The LO vacuum alignment is obtained by 


( 2 ) 

minimizing e<,a!c, 


whose solution is given by po = po with 


, 0 ) 


dpo 


= 0, 


po = arctan ( —etan - ) . 

With this value of p o, the LO vacuum energy density, normalized to 0 at 9 = 0, is [5] 


(A.6) 


(A.7) 


p (2) = F 2 


M 2 (0) — M 2 {6) 


(A.8) 


where M 2 ( 0) = 2 Bm, and M 2 (6) is given in Eq. (10). The NLO vacuum energy density is 


411 = 4ac loop) - M\0) \l 3 + l 7 


(1 — e 2 ) tan(d/2) 


_ 1 + e 2 tan 2 (6/2) 
The perturbation p\ due to the NLO terms is then determined by 


de 


(4) 


dpo 


= 0 . 


(A.9) 


(A.10) 


From Eq. (A.6), it is easy to see that the l 3 term does not have any effect. The vacuum alignment 
is equivalent to removing the tadpole of the neutral pion which causes vacuum instability [13] (see 
also, e.g., Refs. [26, 27, 28, 29]). In fact, with tp 0 = (p 0 , the SU(2) LO chiral Lagrangian does not 
have any term with odd number of pions because such a term is always proportional to 


9 . .9 

cos — sin ipo + e sin — cos <po oc 


de {2) 

O' Cyac 

dp 0 


(A.11) 


This implies that the one-loop diagrams for producing a 7r° from the vacuum as shown in Fig. 3 
have a vanishing amplitude. Thus, we have (d 4a C lo ° P ^^/d po ) = 0. This can be checked 

explicitly with Eq. (A.3) noticing that M 2 (9) = M 2 \^__^. Therefore, Eq. (A.10) leads to a 
solution p\ = pi with 

(A.12) 


Pi = 


F 2 


[l + e 2 tan 2 (0/2)] 
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This is the NLO perturbation to the LO vacuum alignment. Expanding it around 6 = 0, we 
reproduce the result derived in Ref. [29] 

= + (A.13) 

This perturbation produces a CP violating three-pion vertex [28, 29] by substituting Uq with 
if = (po + (pi into the LO Lagrangian, which turns out to be of O (p 4 ). Such a vertex contributes 
to the Trn scattering from O (p 6 ) and to the pion mass only starting at O (p 8 )- Furthermore, 
terms with even number of pions are CP conserving and receive contributions with an even power 
of <pi, so that the (pi-induced terms in the Lagrangian also start from O (p 6 )- Therefore, it is 
safe to make the vacuum alignment at LO for our calculation. It is for the same reason that the 
topological susceptibility up to NLO in the chiral expansion calculated in Ref. [25] agrees with 
that in Ref. [30], where the vacuum was aligned by minimizing the vacuum energy at LO and 
NLO, respectively. 
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